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SLOW WAVE VANE STRUCTURE WITH ELLIPTICAL CROSS-SECTION SLOTS, AN ANALYSIS 

Henry G. Kosmahl, Fellow IEEE 1 
Analex Corporation 
3001 Aerospace Parkway 
Brook Park, Ohio 44142 

Abstract — Mathematical analysis of the wave equation in cylinders with elliptical cross-section slots was 
performed. Compared to slow wave structures with rectangular slots higher impedance and lower power dissipation 
losses are evident below a certain value of kh. These features could lead to improved designs of traveling wave 
magnetrons and gigahertz backward- wave oscillators as well as linear traveling wave tubes with relatively shallow slots. 

INTRODUCTION 

The slow wave, slotted vane structure with or without an opposite parallel plane (Fig. 1) has been studied in 
numerous studies (e.g., by Watkins [2] and by Collins [4]) in more depth. In either a linear version, as in Fig. 1, or in 
circularly bent configuration this structure has been successfully used as slow wave circuit in traveling wave crossed field 
amplifiers or oscillators. The structure has a forward wave fundamental with a backward wave as first space harmonic. 
The latter has served as the most successful circuit for very high frequency: 300 to 2000 GHz, milliwatt power (backward 
wave) voltage tunable oscillator. On the positive side, the favorable feature of this structure is its very high beam 
interaction impedance, when the beam hole is placed into the vane just below y = 0. There the electric field decreases 
slowly as sin k(y - y 0 ) and sin ky 0 ~ 0, y 0 being the center of the beam hole. For this reason the transverse beam slow 
wave structure — coupling coefficient is very large, A/ 2 ^ 0.9 , much larger than in helices, coupled cavities, and similar 
structures with exponential decay. On the negative side, this low pass filter circuit has a typical kh = a Me versus fi 0 L 
characteristic: the curve rises rapidly with high velocity Uph * c/2 to a point where /? 0 L is approximately 7t/3, bends over 
into an almost horizontal (parallel to (i 0 L) shape where the group velocity u g is a small fraction of c, u g ~ c/100. One 
is thus forced to work inside the "bend” region and compromise between high beam coupling impedance K^\/u g and 
losses, also proportional to 1 fu g . With rectangular slots it is not too difficult to achieve K 555 200 to 300 £2 with acceptable 
losses on a short circuit. High impedance makes this circuit also attractive for low power, low voltage applications at 
30 GHz, providing that construction of a circuit with vane spacing equal to approximately 1/5 mm ~ 200 pm is successful ! 
If it is, electrostatic beam focusing with a converging-diverging beam without external fields is feasible. 


1 Work performed for and supported by the Electron Beam Technology Branch, NASA Lewis Research Center, Cleveland, Ohio. 
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The analysis described by Collins in [4] does not assume a constant electric field in the gap at y = 0. By expanding 
the electric and magnetic fields across the gap of the form 

E z (z,y = 0)e ia *= e~' P " Z e ,(0f 

n =— oo 

inside the gap 0 < z < 8; E z = 0, 8 < z < L\ an infinite by infinite determinant is obtained whose zero value yield the 
coefficients E n \ analysis performed^ shows that the zero coefficient E m (as assumed by Watkins [2]) alone gives an 
accuracy of more that 98 percent, that is E z = constant to within 1 to 2 percent across the gap! 

In this study the replacement of the rectangular slot by one with an elliptical cross section with the small semi- 
axis = 8/2 and the large = h has been examined with full rigor, except for assuming the constancy of E z across the gap 
as discussed earlier. Results confirm expected behavior: due to smaller volume and shorter boundaries the slot energy 
and the power dissipation in elliptical slots are smaller than those in rectangular slots. That means doubling the impedance 
and cutting the losses, a very attractive feature for potential applications. Fabrication experiments have shown that 
elliptical slots are not feasible to construct with an aspect ratio h!h = 15, as required by low voltage requirement V 0 < 6 
kV for a 2 W, 30 GHz forward wave TWT. On the other side, higher voltage, low frequency crossed field amplifiers could 
much benefit from higher impedance and lower losses. Since noise jitter and noise figure are proportional to the current 
(at least in the first power), efficiency and noise could be improved. 

The vane structure with elliptical slots could be used for design and construction of a BWO with lower starting 
oscillation current due to higher K and lower losses. 

In the analysis to follow the wave equation of an elliptic cylinder is solved, assuming a single component for E z 
that behaves like the elliptic sine, in analogy to the treatment by Watkins. Although rigorous expressions were derived, 
it would require extensive programming of Mathieu functions to established rigorous results. Such a project should, 
perhaps, be carried out if interest in the possible applications is demonstrated. 

LIST OF SYMBOLS 

A* coefficients of expansion in Mathieu functions 

a , b characteristic numbers in Mathieu equations 

m m 

C constants 

* Unpublished work at LeRC. 
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c 


C em &q), **<§*). Fe y m &q), GeyJ^q) 
ce m (r\,q), se m (r\,q) 
d 
E 


f 

Kq, kh) 

g = (c 0 /V2) ^cosh2^ -cos2r) 

H 

h 

K 

k = co/c 

L 

M 


m 

m 

P 

q = k 2 c 0 2 f4 



u ph’ u g 

W 


W 


X, y,z 


speed of light = (s 0 ji 0 )“ 1/2 
half confocal length of ellipses 
Radial Mathieu Function 

"elliptic cosine" and "elliptic sine" Mathieu Functions 
distance between top of vanes and opposite planar conductor 
electric field, [ Vfcm ] 
frequency 
function, Eq. (33) 

metric factor in elliptical coordinates 

magnetic field, [A/cm] 

height of vanes 

interaction impedance [Q] 

wave number in vacuum 

period in slow wave structures 

longitudinal or transverse beam-wave coupling coefficient 
order of elliptic functions ce m , Ce m , etc., interger 
number of slot in periodic structure, interger 
power [watts] 

parameter in Mathieu equations 
Skin effect surface resistivity [Q] 
phase or group velocity, [cm/sec] 
stored energy, [Joule] 
width 

carthesian coordinates 
^ l t ‘ wave impedance in vacuum [Q] 
ordinary Bessel Functions J n or N n 


(X3 cc 5 ... 

ft. 


GREEK SYMBOLS 

coefficients, Eq. (13) 
slow wave propagation constant, axial 
slow wave propagation constant, radial 
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5 

^n,x 


clear space between vanes 
coordinates of an elliptical cylinder 
£ 0 permittivity of free space 

£(£,r|) = 'F(^) . (J)(r|) product solution of wave equation in elliptical slots 



free space wavelength [cm] 


permeability of free space 


SUBSCRIPTS 

D 

discipitated 

m 

number of slot, interger 

m 

order of elliptic (Mathieu) functions 

n 

summation index for slow waves 

e 

refers to elliptic 

r 

order of coefficients, A r interger 


WAVE EQUATION IN ELLIPTICAL COORDINATES 


In solving the wave equation inside the slots of elliptical cross section we are dealing with the geometry of an 
elliptical cylinder that extends infinitely in the x-direction ( djdx = 0). When the two-dimensional wave equation in 
vacuum 


By 1 


+ 




dz l 


+ k 2 ^ = 0-, * = 


( 1 ) 


is transformed from the cartesian coordinates x,y*z into confocal elliptical coordinates [1], one obtains, see Fig 2, 


3$ 2 3n 2 4 


(cosh 2^ - cos 2t|) c, = 0 


( 2 ) 


where ^(S,, r]) = y(c,) <()(T|)and ^ is a function of £ alone and <J> a function of t) alone. Then we obtain 


<t>L^ + yLf. + — C Q (cosh 2 \ - cos 2 t])\|/-<|> = 0 




3n 


4 


(3) 



and <; = \|f •<)> is either E or H. 

Dividing Eq. (3) by »J/'4> leads to: 


1 d 2 y 
T + 


k c' y l a ^ ^ 

— 2 - cosh 21 f + — t - 2 - cos 2t] = , 

4 * 4> dn 4 


where a is a separation constant. Rearranging leads to two ordinary equations: 


d^ 

w 


2k 1 c 2 


a - 


—cos 2 ti 


<(» = 0 


d 2 y 


2it 2 c 2 

—cosh 2q 


y = 0 


2 


(4) 

(5) 

( 6 ) 


Equations (5) and (6) are called the canonical forms of Mathieu equations. If in Eq. (5) we write ±i£, for T), it is 
transformed into Eq. (6), while (6) is transformed into Eq. (5) if ±rq is used for £. Sometimes this is considered a fluke— 
but a lucky and useful one. Two functions which are solutions of Eqs. (5) and (6), respectively, for the same values of 
a and q = k 2 c 2 /a form a product that yields the desired functions q = if/ # - H ox E. Since a may have any value, 
the number of solutions could become unrestricted. However, the electromagnetic field is a periodic and single valued 
function of t| and £ and its solutions are linear functions. This is possible only for discrete values of the separation 
constants a , called “characteristic numbers” a = a m , b m . They yield, in turn, ordinary and modified Mathieu functions 
of integer order m, corresponding to a m , b m . Observe that we are interested only in cases where q = k c 2 /4 is positive, 
q > 0. 

The periodic solutions of Eq. (5) of first kind are denoted either ce m {r],q) or se m (i],q)—an abbreviation of 
“cosine-elliptic” or “sine-elliptic” — are the ordinary Mathieu functions, while the “Radial” solutions 
Ce m (^,q) and Se m (£,q) and Fey m (£,<?) are called modified Mathieu functions of the order m. They have to be 
combined to form the product solutions that belong to the separation constants a m , 

CeJZ,q)ce m (T),q), Fey m (Z,q)ce m (r},q) (7) 

and 

Se m {^q)se m {^q), Gey m {$,q)se m (r],q) ($ 

that belong to the separation constant b m . 
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The corresponding functions in cases of circular cylinders (e.g., helix) are the ordinary or modified Bessel 
functions (e.g., / m (v m x)^(7 m ^) or J m {y m ^)N{y m x) and cos(mri) or sin(wn). Notice that, in general, a m * b m and that 
mixing of Mathieu functions other than that shown in (7) and (8) is not permitted. The power series for ordinary and 
modified Mathieu functions were computed and tabulations may be found in [1] and [3]. 


DISPERSION EQUATION 


We shall derive in this paragraph the dispersion relation /J(oo) = f{co/c) = f(k) for the finned structure shown 
in Fig. 1 with slots having an elliptical shape. Since it is assumed that the slots extend infinitely in the x-direction 
(perpendicular to the plane of paper) we are dealing with elliptical cylinders with no x variation, d/dx = 0 ; it is very 
helpful and necessary to invoke the comparison with the case of rectangular slots, as treated by Watkins [2]. Watkins 
takes only the O-order case with ns variation of E and H along the boundary y = 0. The appropriate solution for standing 
waves in cartesian coordinates is then [2] in the m-th slot (not to be confused with order in above) 


£ sin fc(y+ft) e ~i$ 0 mL 
z ° sin kh 


mL < z < mL + 5 


(9) 


E = O, mL +8 < z < (m + l)L 


„ . I £ n cos k(y+h) _,n m L 

= - curl E\H = i M- r 0 

x smkh 


(9 a) 


( 10 ) 


H x = 0, mL + 8 < z (m + l)L (10a) 

e,-b,.h j -h,.0 01 ) 

We take notice that both E and H are independent of z inside the rectangular slots! We are dealing with a TEM wave only . 
Returning now to the case to be treated in elliptical cylinder coordinates it seems appropriate to use the Ansatz yielding 
results of the form of Eqs. (9) and (10). Since H x is in the direction of the axis of the elliptic cylinder that has a maximum 
at the apex (bottom of slot), r| = 0, we put the linear combination 


H x(^) = Cce \ (O’tf)' 


+ Eey^q) 


( 12 ) 
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Equation (12) is the (product) solution of wave Eq. (2) in elliptic cylinder coordinates. Note that |) is independent 

of*. cej(T|, q), Cgj (^,9) and Fey\ are the angular and radial solutions of Eqs. (5) and (6), respectively and are given in 


[1] and [3], 


= cost] -- 


l + £ + £- 
8 192 


+ -• 


cos3r( + ^2^ + £ + •••jcos5r| - 


(13) 


For values of q not too large (q < 8), the series con verges quickly. In our case q < 7T 2 /l 6 < 0.61685 The radial functions 
are 




7t 

l’ q 


M 1 


-/tlyj^T^cosh^j + cosh - A ] 5 J 5(2 ~Jq coshZ^ + - 


(14a) 


ce. 


Fey^q) = 


-AfjV, (l-Jq cosh + A3 N 3 {2^[q cosh 2;) - Aj N 5 ( 2 ~Jq cosh (;) + (14b) 


In Eq. (14) 


ce. 


fn ") dce x {r\,q) 


\+-q 

( 2 \ 

i + £ + -£- 

+ -L4 I + a + 

W q ) dn 

ti=n/2 

8* 

l 8 192 J 

192 t 6 ) 


After dividing Eq. (14) by a] , Eqs. (14) may be written 

J * 

cc 

Ce x {%,q) = - 


Ce '\J' q 




A 1 A 1 

( 2 -sfq cosh jJ^Jq cosh^) + —jJ 5 ( 2 -yfq cosh - 


(13a) 


(14 a) 


A* 


Fey { {%,q) = ~ 


ce \\2' q 




N^yfq cosh^j cosh (;) + -. . . 

A l 


For an ellipse whose long half axis a — c 0 cosh^ » b — c Q sinh c, - — — half gap width . 


= Ja 2 -6 2 = Ah 2 , cosh£, = — = 

" 4 c„ 


. 2 -«i 


hi 

Lii 


4h 2 


1 8 2 
~ 1 + 7T 

8 h 2 


(14b) 


With - = -L as required for a 30 GHz TWT, cosh| 0 = 1.000556; and £, 0 = 0.03345, 
h 15 

with a small error cosh £ can be taken as 1 in the arguments of the Bessel functions. Thus: 
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(14c) 


Ce x {%,q) = - 



^f y 5 K) - 


And because S « h 7 c Q ~ h, and Eq. (J4b) may be changed to 


ce , 


Cefaq)* — 


K 

T q 


Fi 


J \ ( kh ) --T' J 3( kh ) + -T’ J 5 ( kh ) 


and similar for the Fey\ function. 


(14d) 


Expressions for the coefficients Ajn+x derived later. On the long axis £=0, on the ellipse % = and the argument 

of the Bessel series, kh> changes to 1.00055 kh with negligible change in the Ce and Fey value. Thus, is almost 

independent of £, or correspondingly, HjJy) is independent of z in rectangular slots (Eq. (10)). The expression for the 
electric field components ^ an ^ follow from Maxwell equation curl H = +icae o E 




kg 


Ce l {^q) + ^Fey ] (^,q) 


= ifclit = 


dCe,(^q) | C 2 Fey,(^q) 




* 


(15) 

(16) 


with g = (c 0 j -J2 j ^/cosh 2^ - cos 2rj , the elliptical metric factor. The arc lengths ds\ and ds 2 , the hyperbolic and elliptic 
arc lengths, respectively are 


ds x = gdt; and ds 2 - gdr] (17) 

Note that ds\ being perpendicular to the elliptic contour £ = const corresponds to changes across the gap (dz in [2]) and 
ds 2 is perpendicular to hyperbolic lines t] = const (dy in [2]). Differentiating Cei(£,#) with respect to £ yields from 
Eqs. (14a) and (b): 


dCe&q) 

X, 



■ 2-yJqsinha 


9y,(p) A* dJ 3 (p) 

dp Aj 1 dp 


+ ••• 


(18a) 
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dfoyi (£,<?) = 


n 

ce i\ T><? 




* 2^sinhd 


8^,(p) Aj dj^jfi) t 

dp A, 1 dp 


(18b) 


with 


P - 2V?cosh«, d -iM = /» - = /,w - ^,(p) 


. However, the tangential component of E is Zw, on 


etc. Because, in the 30 GHz case, sinh£ = 0.03345, | | « 

the elliptical surface Eq. ( 1 8) shows that = 0 only for £ = 0, that is on the confocal line. Because £ « 0, sinh£ ~ £ 
is very small but not zero. To force E r | to be rigorously zero on we have to force the brackett of Eq. (16) to become 
zero on the surface £ = Thus: 


* 


+ 


S 0 


C 2 dFeyj^q) 

C, 31; 


= 0 


So 


(19) 


Note that when £ = 0, the argument of N becomes simply 2 yjq = kcjl , regular and finite. 

The constant C( 2 ) is determined from the requirement (Eq. (19)) 

Observe that the Bessel expansion in Eq. (18b) is the same as that for Ce \ given in Eq. (14a) since both, the J n and N n 
functions have the same recurrance relations. It then follows from Eq. (19) 


or 


V 3 


'(i) 


A 

J \( kc 0 cosh ^) - ~A J l{ kC o C0S H) - +• 


+ C 


( 2 ) 


A 1 A 1 

JV^cosh^) - -j-W 3 (k 0 cosh^) + cosh^) -+■•• 

A i A i 


= 0 




( 20 ) 


fe^sinh^, C (1) 


A 1 A 1 

J[{kc 0 coshE, 0 ) - - 4 - [kc 0 cosh ) + -f J' 5 (kc 0 cosh^ 0 )~ + 


+ C, 


( 2 ) 


A 1 A 1 

N{{kc 0 cos^) - ^(fcc 0 cosh^) + -4-N' 5 (kc 0 cosh -+ 


= 0 


( 21 ) 
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In Eq. (2 1 ) the prime over the Bessel functions designates the derivative. Equation ( 2 1 ) is not an eigenvalue problem 

and it must be satisfied for all frequencies, that is regardless of the value of kc 0 = (a >fc)c 0 . Since, as discussed earlier, 
sin h^ 0 * 0 , the expression in the wave brackets must vanish. 

We have then, using p Q = kc Q cosh^ = kc Q 


C 2 (q) = -C x (q) 



The coefficients A[, A 3 , ^A' 2r ^ x may be obtained from recurrence relations as given in [1] and [3]: 

The relation for the desired function ce l ( Tj,q ) and r > 1 is 


( 22 ) 


(flj - 1 - q)A\ - qA\=0 


(23a) 


a 


l 




(23 b) 


For q < 1, the series converges rapidly and the required ratios A^ r+ 1 / Aj 1 are easily obtained for a given value a \ — the 
characteristic number (separation constant) which are listed in [1] and [3]. The required expression for a\ is 


Then 


1 


2 3 

i h <r 

= 1 + <7 - 

H 8 64 1536 



(24) 


(25) 


Table I lists the values of a x , a] , A^ for a range of kc 0 values as parameters. As q changes with frequency, the ratios are best 
computed as functions of q. Nevertheless, for ^ < 7 t 2 /l 6 , A-^/aJ = qj 8 ; A^j A* = q 2 /l 92 to demonstrate the rapid 
decay of coefficients A 2r +i with increasing r. The derivatives of the /(p) and N(p) functions with regard to the 
argument are([z(p) either J(p) or N(p)]): 


^ Z n(P) _ _SL 

dp p 


z » + z *-i(p) 


(26) 
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Equation (22) may then be evaluated. As an example for kc Q ~ kh = 1.4 and = 0.03345, p = 1.4, 
one gets C(2) = - 0. 1 85 1 C(i). One should keep in mind, however, that C(2) = C(2>(<y), that is C(2) depends on frequency. 
This is also true of the coefficients A. r {q) — A r [k 1 / 4 j . Thus, the C(i), C(2), and the A r ’s have to be computed as function 

of q. Table II shows the dependence on q of some important parameters and Bessel and Neumann functions to 
demonstrate their behavior. Note that for q < 1, that is of interest in this study, |a^| « |a^| « |aJ | and the Bessel- 
Neumann series in the Eqs. (20), (21), (23), and others converge rapidly. Nevertheless, the expressions are cumbersome 
and for practical evaluation programming becomes indispensable. 

We are turning now to the important expression for ^(^77) that has to be matched at the gap y = 0“ or 7] - kh 
to the E z component of y = 0+ From Eqs. (12) and (15) 


7 7 


Ce^q) + Feyfcq) 


"(l) 


(27) 


and cel(ri,q) = dcc } (77, q)jdr\. We have to normalize the parameter C(i) = Qi)(<?) such that for T]-kh: 

-i^ce[{kh,q) c m [ ] = E g 


(28) 


The [ ] bracket designates the Bessel functions summations in Eq. (27) 


kc n 


■"(l) 


= 1 


Z 0 ce[(kKq)\ } 


(28 a) 


At the center of the gap £ = 0. The value of 77 = kh at y = 0 is close to kh = nj 2 . Let us put 

77 = kh = — - 0, with 6 « 1 . Then : 

' 2 


g — ^/cosh 2£, — cos 2r| 



and kg in Eq. (28) becomes kg = kc 0 . As discussed in the introduction the assumption E - E 0 in the gap is better than 98 
percent accurate for narrow gaps and h » 5. Using C(i) from Eqs. (28) and (12) one gets 


H x &ri) = -i 



(^hq) 

ce[{Tl,q) 


(12 a) 
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1-1*1 
8 1 


( 


-cef(ri, q) n=kk = siRkh 

= sin kh\ 


% Q Q 
l + .z.+-z_+. 

8 192 




v 


sin 3 kh 


sin kh 


192 l 6 J s: 


sin 5&A 


sin£/i 


, „ sin 3 kh „ sin 5&/i 

1 - 3a. + 5a, +• 


sin kh 


sin kh 


( 29 ) 


Thedispersionrelation(Eq.(1.67)or(1.79)in[2])isobtainedinourcasebymatching H+ to H x asgivenby Eq. (12a) 
at the boundary of the gap, given by T] = kh . Using Eq. (1.65) [2] 


* 7 0 Zw v 

^ n=-oo *n 

and taking, as discussed earlier z = 5/2, the center of the gap (£ = 0) and H~ from Eq. (12), one gets: 


(30) 


ce^kh,q) _ y A^ n 8 


(cot kh)f(kh , *) = | C -- ^ - 
1 sin kh 


, cos3fc/i 


i COS 5ta 


l-a,(*) +a<(9) 

3 ^ cos£A 5 ' cos kh 




5 ^ 


, „ sin 34A ^ sin 5#* 

1 - 3a. — + 5a, + 

j hi> 5 - 


= £_ y 

L y c 

J\ Oa 0 


(31) 


(32) 


sin £/i sin kh 

where a 3 (#), a 5 (#)'•■ are the coefficients of the cosine expansions given in Eq. (13). We now rewrite 
Eq. (32) 


5_M^ 

j + y M n! M o 

m 

[j , , M xl M n , M_ 2 /M 0 ( M 2 /M n I 1 

^ Y oe 

n^ ne ^ oe . 

L Yoe 

7_,/7 0 Y-i/Yo Y-l/Yo Yl/Yo 


For large n, M n = (sinn7r*5/L)/(rt/r*5/L) — > 0 and y n = [5 n ItdiS/L — > °o 
Further rearranging of Eq. (33) gives 


y 0 e c o~y 0 e h = -[ M oe 


MnJ M a 


l+y ■ 

“o Inehoe 


tan kh 
f{kh,q) 


(33) 


(34) 


The subscript "e" denotes quantities belonging to "elliptic" solutions. The function f{kh,q ) is equal to 
1 as k = CO/c -» 0 because q = k 2 c 2 /4 goes to zero as ca 2 and a 3 , a 5 , ■■■ also become zero. It also equals 1 as 
kh — » k/2 and has as shallow maximum around 1.08 at kh-nj 4. Around the selected design point 
kh ~ 1.45, / = 1.01 ~ 1 . Comparing Eq. (34) with Eq. (1.67) [2] 
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_J_ = k . — 

tan kh L ^ y n L y „h 


s u m - 

V 


or 


»‘T 


kh tan kh M, 


1 + y KdMo. 

h tn/y 0 


MJMo 


i + Y- 

r »/ro 


(1.67) 


(1.67) 


and assuming that the sums in the square bracket are not much different from each other in both cases one obtains an 
approximate expression 


/ 

1 

r — J 

L 2 - — 

O y 

4 

\ 

/ 


y o e=Vo 


1 


From Eq. (35) it follows 


h L - — 


g 2 f(kk,q) ■ kh 


(35) 


^oe = , * + 


2 . ft 2 rl 

c 2 0 f 2 m 2 


(36) 


The plots P Q Ljn and p oe ■ L/k versus kh as ordinate are shown in Fig. 3. Expression (36) is a good 


approximation when rj = kh = ;r/2,then cos 2 t| = cos it = — 1 and g = c 0 /V2 -y/cosh 2^ + 1 = c g cosh c, = c 0 .The 
case kh - njl is of interest for low voltage designs because then h = A 0 /4. When ti is much less than Jt/2, then 


g = -^cosh 2<^ - cos 2 fcA = —j=- — cos 2kh = —j= -^2 sin 2 kh = c 0 sin kh 

because 0 < ^ < « 1 and cosh 2^ 0 » 1 + 2 ^ * 1. In the general case, the matching of H~ to H+ at £ = 0 is 

accomplished by the expression 


, cos 3 kh cos 5kh 

_ 1 - a, + a s +•• 

, , 7 -x. cosWz cosfc/i 

i — “ (cot khjkg * 


sin 3&/z _ sin 5kh 

1 - 3a, + 5a< 


E 8 M n 
= i — k— > — 2- 
Z L ^ Y n 


sin A:/z 


sin/:/i 


or 


2 


tankh 


/ \ 5 'V"’ 5 A/ V ^ 

***•«)- 

alln '« n*OV 


M./M„ 

! H £L/ — 2. 


At £ = 0, g = -^=--^2 sin 2 kh = c 0 sinfc/i, Eq. (37a) becomes 


(37) 


(37 a) 
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6 Af f A/ / Af rt 

/(w,? =-— y i + , £ 

sin ^ L y 0 £k{ yjy 0 


c 0 sin kh - cos kh 


(38) 


Equation (38) is always accurate when ^ = 0 at the center of gap. As mentioned earlier, a rigorous evaluation of 
Eq. (1.67) as described by Collins [4] indicates that the M 0 = [sin P 0 (5/2)]/[P 0 (5/2)] should be replaced by which, 
for small gaps, amounts to 2 percent correction. Thus, the corrected formulae should be rewritten as 


1 

mkh if y, t „/ y „ 


kh tan 


(1.67 a) 


for rectangular slots [2] and 


c a cos kh - 


1 M 5y^ 

/(“■«) y« 't-tis 


1 + 


'Yne/'Yoe J 


(38 a) 


as more correct expressions for the dispersion relation. 


Power Dissipation in the Elliptical Slot 


The expression for the power dissipated on the walls of the elliptical cylinder is given by 
P D = gdr\dx = ^ Sdr\ 


(39) 


[^[£2] = 34^j p/X\ p[£2cm], X.[cm]) : 

For small values « 1, the ratio Ce l (^ o ,q)/Ce x (0,q) in Eq. (12a) may be taken as 1 with an error of 
0.2 percent in the worst case and 


"xfod = 


frY 


\ Z o J 


1 o) [ce{(kh,q)f 


(40) 


With g = cJ4 2 -JcoshlZ, - cos2n Eq. (39) becomes at £ = % 0 


p -4 

P D- 2 


’ f T- ^ 2 


Jo, 


(kc 0 f l\ce\(kh,qrf ■ c'fa.qyfo osf.2^ - cos2, dq (41) 
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(41a) 


2 


|XY 

'( K } 

2 wc o f 

Z *J 

[ce{{kh,q)j 

V2 J( 


The expression - cos2ri/cosh2^ 0 leads to elliptic integrals which are well tabulated or easily computed, but do not 
yield a useful closed form formula. To examine the behavior of the integral in Eq. (41) let us rewrite 


^cos/i2£ 0 


cos 2 77 = Jl + 2%q + -Z 0 - 1 + 2 sin 


77 + ... 


Now the ratio — , the minor to major half axes, equals to 
a 


b ^512 
a h 


sinftg 0 

cosA^ 0 


= tan/i^ 0 . 


(42) 


For the 30 GHz Forward Wave TWT at 6kV: 5 / 2h = 1 / 30, and = 0.0333. For a high power Forward TWT 

magnetron amplifier < 0.06. The highest ratio 8/2/t ~ 0.25 appears in the designs for 1THZ BWO, with the 
resulting value = 0.25. The value of the intergral in Eq. (41) is then approximately 10% larger than that obtained 
with negligible % 0 . Due to the smallness of ^ the integral and the integrand in Eq. (41) will be different only when 


sin t| < %q. Figures 4(a) and (b) show computed values of F(x,^) = J Q ^0 + s * n2 “H dT| with £ 0 =0.25 and 
r x { x 

0.0666 and I(jc) = sin r\ dr| = - cos v\ . 

Jo *0 


For £ 0 = 0.06, at x = 0.2; F(j t,£) = 0.025062 and I(x) = 0.01993342 or an error of 20 percent. But at kh = 1.5, the 
actual upper integration limit, F(x,^) - 0.94 anerrorofonly 1 percent. Thus, in computing the total loss we are justified 
to neglect in Eq. (33). Then: 


1 E } rxA ? 

p D = r k -f wc 0 L cef{T],q)smr]dr] 

2 \ Z o) J ° 


(43 a) 


Using Eq. (13) for q mm ^ 0.5625; q 1 2 * < 0.316; etc. we may neglect higher powers in q and a. ce^(rj, q) then 


becomes 


ce ! 2 ( r l><?) = cos 2 ri 


1 - a. 


cos 3t| cos 5r| 


n2 


COS T| 


+ a 


- + -* 


COS T| 


= COS T| 


1 - 2a 3 ^4cos 2 r| - 3 j + 2a 5 (cos 4 T| - lOsin 2 t^cos 2 J] + 5sin 2 r| j + — j 

= cos 2 r|^l - 8a 3 cos 2 T| + 6a 3 j 


(44) 


inserting Eq. (44) into the integral in Eq. (43a): 
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( 45 ) 


J cef(r\,q) sint)^ = cos 2 r] sin r|(l - 8a 3 cos 2 r) + 6a 3 )dr| 

^•(l + 6a 3 )|*Q + |a 3 cos 5 ti|*q = 1 + ^ a 3 ^ _ cos 3 + |a 3 (cos 5 kh - l) 


cos 


Integration from 0 to kh gives the approximate power loss on one side. Using Eq. (43a) for two sides one gets: 


P De=T 2 ^ 


-5l 

3 


' f £ \ 


2r 


L 




K Z o, 


kc„ 


V 


kc„ 


K ce[(kh,q) 
-|2 


wc. 


1 + 6ou 


|l - cos 3 fc/ij + — a 3 ^cos 5 kh - lj 


ce. 


[{kKq) 


wc 


(l + 6a 3 ^l - cos 3 kh} + -y^a 3 |cos 5 kh - ijj 


(46) 


From Eq. (29) 


ce[(kh,q)Y ~ sin 2 kh 


* _ sin3fc/t . sin 5kh 

1 - 3a, — + 5a< +■ 


l2 


sin kh 


sin kh 


•y 

sin kh 


, , sin 3kh 2 sin 2 3kh^ 

1 - 6a, , ,, + 9a 3 — = — 
sm kh sin kh j 


(47) 


if a 5 « a 3 as well as products a 3 a 5 are neglected. Now, at the operational point of interest 
kh =1.4 to 1.5 and sinUi = 1; sin 2kh ~ 0.1 and sin 3kh = -1 ; thus, for kh = 1.5, Eq. (47) becomes 

[ce[(l.5,^)] 2 = (l + 6a 3 + 9a 2 ) sin 2 kh 


and Eq. (46) becomes approximately 


P 1=3 

r De I 




K ) 2 


2 

sin kh 




(48) 


The corresponding result for rectangular slots of equal width w and height h and gap 8 is 


P DR 


R h 


sin 2 kh 



wh\ 


1 + 


sin 2 kh 6 

2 kh 


Equating Eqs. (48) and (49): 


n2 kh 5 
2k h~ + h_ 

j( KC oJ wn ^ 1 4h 2 - 2 sin 2 kh (l -4.8a 3 ) 



(49) 


( 50 ) 


or 
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kc„ 


l 


(1 + 


sin 2 kh 5 
2k h +h 


^(l-4.8cc 3 ) 


= 1.53 


(50 a) 


For the ratio P DR jP De it then follows from Eqs. (46) and (49) 

5 sin2tt i - 6cx 3 + 9<x | 

Pqr m 3 h h 2 kh sin ^ s™ kh 

P De 2 c o (^ c 0 ) 2 (l + 6a 3 )|l - cos 3 kh ) + ^a 3 (cos 2 kh - lj 


Thus, below approximately kco~ 1 .53 the losses of an elliptical cylinder are smaller than those of an equivalent 
rectangular counterpart. Moreover, they decrease roughly as [kc^ = k 2 h 2 \\ - (5 2 /4/z 2 )j decreases. For a given 
frequency k it means either a decrease in h (higher voltage) or increase in 8 (gap width). The former is of interest to 
traveling wave magnetrons, since reduced losses result in lower noise output and for BWO’s at harmonic frequencies, 
where increasing hlh is a necessity. 

Since the a coefficients are functions of frequency, accurate evaluation of Eq. (46) and others requires extensive 
programming and computation work. The power losses in the space above y = 0 are identical to those obtained in 
rectangular coordinates. 


Stored Energy in the Elliptical Slot 


Since H x {^,rj) is the only magnetic field component, the stored energy is given by 

-w; Jo” jn*« 2 *~ <52> 

because integration from 0~4o covers only one half of the elliptical slot. 

For small £ « 1, we neglect the dependence of Hon £ as discussed earlier (see Eq. (14)) and obtain 




2 kh , J 

J J £ M 2 


4=o-n=o 


J 


— ^’^ -(cosh 2S, - cos2r|)oiqrfri 
\ce[{kh,q\ 


(53) 
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/ \ 2 ✓ \ 2 2 - 

= ^- — - —2. — w sinh 1; 0 cosh 1; 0 f jcos 2 T|(l + 6a 3 )-8a 3 cos 4 T| + 

2 K Z o J K ce i ) 2 L J ° 

2a 3 (cos 4 ri - lOsin 2 ricos 2 r| + 5sin 2 t|)J«fr| - cel(r\,q)cos2r\dx\ (55a) 

In developing Eq. (55a) higher powers of 03 and as were neglected for simplicity, 
a 3 being < 1/16, a 2 < 1/256 and a 5 < 1/192 in the worst case. When integration of the terms in Eq. (55a) is 
performed, the following results: 



Except for having neglected higher order terms of a, Eq. (56) is an accurate expression for W e . 

Recall that c g sinh £ 0 = 5/2 = half of the minor axis and c g cosh^ = h = half of the major axis and that 


2 

f^° f f ^-(cosh2^ - cos2r|)d^dri 

Jo Jo Jo Jo 2 

2 ^0 2 

= — J cosh cos2t| dr\ (57) 

0 0 

c 2 

= JL. s\nh^ 0 cosh^ 0 * 2 n = 7T-/i-5/2 = area of an ellipse, because the second integral in Eq. (57) is zero. Let us 
examine the result in Eq. (56). Since sinh and cosh may be expanded in a series and for £ « 1 the series converges 
rapidly, Eq. (56) may be presented in the form: 
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W-b. 

e 2 


V Z oJ 


kh 


KL wc 2 £ - 
(“,f ° 5 ’ « 


sin2AA sin4£A* 

1 + 2a, + 

3 kh 4kh J 


( 58 ) 


Let us compare Eq. (58) with the simplest form of stored energy in rectangular slot of height A and width 8 


W, 


= to. 


( T? ^ 2 


v z »; 


5h w f sir^AA'j 
sin 2 kh v 2 kh J 


(59) 


Taking the ratio of Eqs. (49) and (48) one gets approximately 


2 ( sin 2 kh\ 

Wr = 2 H(AA,g)] SA [ 2kh J 


W e (kc 0 ) (kh) c 0 l o sin 

sin 2£/i 

^ I + bcx^ + 

» 2 - 


K>) C A 1 + 2a 3 + 


1 + 6ou + * 

3 2kh 
sin 4kh 


4kh 


(60) 


(60 a) 


But 5/2 = c g sinh % 0 ~ c 0 ^ 0 for £, Q «\, and then 


W D 


W. 


'e {kc 0 ) l + 2a 3 + 


, sin 2kh , 

1 + + 6a, 

2kh 3 
sin 4 kh 


4 kh 


(60 b) 


For kc 0 2 S 2, the numerical value of Eq. (60b) is larger than 1, a significant fact, that impacts the interaction impedance 
toward higher values in elliptical slots. When the slots become less deep, kc 0 = 1, Wr/W c ^ 4 indicating large 
improvements in impedance. The stored energy above the slots 0 < y < d is identical to that obtained in rectangular 
coordinates, designated W u 


n =— oo * n 


sinh 2 y n d 

1+ 

sinh 2 y n d 


(61) 


The total stored energy required for computing of the interaction impedance is the sum of Eqs. (56) and (61 ). Note that 



DISCUSSION OF RESULTS 


The analysis developed in this study took as an approximation only one wave downwards into the half ellip- 
tical slot represented by the ce x {T[,q) function and its derivative cej (T|, q ) . They correspond one by one to the simple 

TEM wave representation {H ~ cos ky\ E ~ sin ky) in the rectangular slot which leads to E z = Eq = constant across the 
gap and to an error of less than 2 percent, as discussed in the introduction. To include higher TM waves in elliptical 
coordinates would require taking all ce m {r\,q) with corresponding Ce m (& y< j) and Fey m (^, q ) functions and obtain then an 
infinite by infinite secular determinant equal to zero. The resulting improvement in accuracy is probably in the same 
range of 1 to 2 percent. Most numerical results quoted in this study were obtained by neglecting €£5 and higher powers 
of CX 3 ^ 3 , OC 3 , . . .j . For q ~ 0.5 the errors are q 2 /6 4 * 1/256 for each term. Altogether all the errors are likely to amount 
to 5 to 10 percent, an error still acceptable for comparative evaluations. 

The virtual constancy of E z across the gap S is evident from 

When h » 5 and 5 
d E z I dz = 0 and E z ~ 


„ dE dE 
V£ = — ^ + — -2- = 0 
oy dz 

« X 0 , E y being tangential to the slot wall is nearly zero inside. Then 


constant. 
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APPENDIX I 


TRANSVERSE BEAM-WAVE COUPLING COEFFICIENT OF THE SLOTTED 
SLOW WAVE VANE STRUCTURE 


Figure 5 shows a cross-section of a vane with height h and width w and a beam hole of radius b « w. 
The center of the hole is located at O' or a distance -y 0 below the top of vanes, y = 0. In the simple (but very accu- 
rate) single TEM wave theory the electric field in the slots behaves as equation (9) 


£ z (y) = 


sin k(h + y) -ip a mL 


sin kh 


( 9 ) 


We introduce now circular coordinates inside the beam hole with radius b centered at y = -y 0 . Consider a 
narrow strip of thickness dy located at y = -y 0 + b cos a . The half width BA = b sin a. The area of the strip dA is 


dA = Ibiinady = -2b 2 sin 2 ada 


(Al) 


cos a = 


Z12k ;sina =l,Jr±iLf 
b j l 1 J 


The area of the hole is, of course 


A = - 2 b 2 f sin 2 a da = nb 2 


(A 2) 


The magnitude of E z at y = -y 0 + b cos a is: 


_ Si. 

1 zW| 0 sin kh 


(A3) 


The transverse beam coupling coefficient M 2 is defined: 


MX 


7zb 2 E 2 ^ 


E 2 dA = 


- 2b 2 


7tb~ sin kh * 


r-y a +b 


Ly o -b 

l b J _ 


sin 2 k{h + y)dy 


2 

/rsin kh 


•C sin 2 ® sin 2 k{h -y a + b cos ajda 


(A 4) 


with 


sinfc^/i - y Q ) + b cos a = sin*(/z - j 0 )cos(/i> cos a) + cos k(Ix - ^)sink(fecosa) 


and, because kb = (kh)— = 0.15 < 1 , we may expand 
h 
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APPENDIX n 

A BRIEF LOOK AT THE SUBMILLIMETER BACKWARD WAVE OSCILLATOR 
IN ELLIPTICAL CROSS-SECTION SLOTS 


To date the only successful B WO in the frequency range 500 to 2000 GHz is the vane structure (analyzed in this 
paper) as used in the first space harmonic with opposite group and phase velocities (see Fig. 6). 

Let us examine the design parameters of such a backward mode for the rectangular slots. The favorable range 
lies at a phase shift pL ~ 2 k - nj 3 because the group velocity becomes much smaller than c and with it the coupling 
impedance increases to required values to permit the start of oscillations. 

The phase velocity v, = co/ P_ 1 is obtained from 


with fi 0 L chosen to be — it follows 


P-l ■ L = 2k — P 0 L 


H ~ l 3 


go coL 



(A 7) 


(A 8) 
(A 9) 


Assume the choice of the beam voltage V 0 = 5000 V. Then the electron velocity u a = 4. 166. 109 cm/sec. For a frequency 
of oscillations / = 1000 GHz = 1012 Hz one gets 


c oL 



5 

3 g 5 _ 5 4,166-lQ 9 

~ U °2 Jtf 6/ 6 MO 12 


3.47- 10 -3 cm 


(A 10) 


The operational point on the kh = F(/}L) dispersion curve has a value kh ~ 1.4 =(co/c)h and 


h = 1.4-^- = 6.684 10 3 cm 
2 nf 


(All) 


Thus h!L = 1.926 and because 8 = L/2 


h_ = h_ 
L 28 



If the slots were designed to have an elliptical cross section with the depth h (half major axis) and c g 


(A 12) 

* JO* - <5 2 /4 = h 


we obtain for the important parameter kc g - kh = 1 .4 and from Eq. (60b) 
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Thus the stored energy of elliptic slots is almost two times smaller, or the interaction impedance approximately two times 
higher. To examine the impact of this result on the performance of the BWO we refer to Figs. 11.1-3 and 1 1. 1 -4 in [5] 
reprinted as Fig. 7. In Fig. 7 the parameter Q/N, the quotient of the space charge parameter Q and the number N of 
electronic wavelength on the structure is shown. Note that in the linear theory Q is independent of I 0 \ 



f CO V 2 
<? 

\ 2co J 


1 

4 C 3 


R 


2 a l_ 




= R A 



4j2 ji 3 b 2 e o Kf 2 


(A 14) 


where C is Pierce’s gain parameter, u) p the unreduced plasma frequency, 01 q = R ' a) p, R being the plasma frequency 
reduction parameter, b the beam radius, and K the effective interaction impedance. Knowing Q from Eq. (A 14) and 
assuming a range of values of N the value of CN = (I 0 * K.\f4V 0 )m ’ N necessary to start oscillations has been plotted in 
Fig. 11.1-3. Knowing C, the necessary current I 0 may be calculated. 

An examination of the curves in Fig. 1 1 . 1 -3 indicates a double benefit for the starting current as result of increased 
interaction impedance K_i : First, increased value for K_i reduces the value of the passive mode parameter Q according 
to Eq. (A 14). For a selected value N and an assumed or measured value of loss decreasing Q moves the corresponding 
Coordinate to lower values. Then, secondly, because = Rco p increasing K-i decreases Io, the starting current, such 
that the product Io K_i remains constant. It appears, therefore, that doubling the impedance would reduce the required 
lo-start by a factor higher than two, a very valuable advantage. As an example, the Thompson-CSF made ITHz BWO 
which utilizes the very same vane circuit with rectangular slots, has the following design parameters at 
f=lTHz= 1000 GHz: 


L = 

40 pm 

8 = 

20 pm 

h = 

75 pm 

N = 

30 

w = 

90 pm 

b = 

22 pm (beam opening radius) 

V 0 = 

7 kV; I 0 = 22 mA; 

Jd 

=1750 A/cm2 (current density in the beam) 


This extremely high requirement for current density imposes great difficulties on the cathode life and gun construction. 
Any relief by factor 2 or more would be highly beneficial and necessary. 



TABLE I. 


kc 

k 2 c l 

« = 4 

a ! 



L4 

0.49 

1.4581 

-0.06508 

0.001356 

1.2 

.36 

1.3431 

-.047025 

.00716 

1,0 

.25 

1.242 

-.03222 

.00034 

.75 

.1406 

1.138 

-.01788 

.00105 

.5 

.0625 

1.062 

-.00787 

.00002 

0 

0 

1 

0 

0 


TABLE n. 





C (2 /C(,) 

Q 

c (l> 




0 

0 

.5 

.2422 

-1.4717 

■ 

.3926 

1.0 


-.7812 


.6272 

1.2 

.4983 

-.6211 


.6433 

1.4 

.5419 

-.4791 


.6305 

1.5 

.558 

-.4323 
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Figure 5. — Evaluation of the transverse beam coupling coefficient 
The 0 of the coordinate system is located at the center top of the 
vane. 
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/-v = c 




Figure 7. — Oscillations are produced in a backward-wave oscillator 
for values of CN equal to or greater than the value give above. Its 
value is a function of the space-charge parameter and the total 
circuit loss. From Proc IRE. 
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